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Abstract 

Anomalies are an infrared effect, but are often realized in effective theories in a non-trivial 
way. We study the chiral anomaly in Soft Collinear Effective Theory (SCET), where the anomaly 
equation has terms contributing at different orders in the power expansion. The chiral anomaly 
equations in SCET are computed up to NNLO in the power counting with external collinear and/or 
ultrasoft gluons. We do this by expanding the QCD anomaly equation, using the tree level (LO 
in Us) relations between QCD and SCET fields. The validity of this correspondence between the 
anomaly equations is confirmed by direct computation of the one-loop diagrams in SCET. 
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I. INTRODUCTION 



In calculations with a hierarchy of scales it is often useful to work with an effective field 
theory, in which ultraviolet degrees of freedom have been integrated out. Examples of this 
are Heavy Quark Effective Theory (HQET) and Soft Collinear Effective Theory (SCET). 
Or as a more well known example, integrating out a heavy quark. 

Although anomalies [ll, Q] come from ultraviolet divergences, they are actually an infrared 
effect because only massless particles contribute Q. If masses are generated through a Higgs 
mechanism as in the Standard Model, fields are massless above the symmetry breaking scale 
and can contribute to anomalies too. 

Since they are an infrared effect one would expect effective field theories to reproduce the 
anomalies of the full theory, but this is not always a trivial matter. As an illustration we 
consider the case of integrating out a heavy quark doublet^. Simply removing it would spoil 
the anomaly cancelation for the gauge symmetries. D'Hoker and Farhi show in a detailed 
analysis that the effective action contains additional Wess-Zumino terms and terms involving 
the Goldstone-Wilczek current, that restore the gauge symmetry jl, HI. 

Another example is the axial anomaly in Chiral Perturbation Theory (ChPT). By itself 
ChPT would not include reactions such as tt" 77 or KK tt+tt^tt^. They do not show 
up because of the symmetry M — ^ — M, where M is the matrix of Goldstone bosons, which 
is not a symmetry of QCD. Again one needs to add Wess-Zumino terms and additional 
terms that couple to the gauge fields f^, f^. 

The goal of this paper is to study how the chiral anomaly is realized in SCET. SCET 
has been introduced to describe the dynamics of energetic light hadrons and provides a 
systematic way to separate the hard, collinear and soft scales [1, 0, [13, It captures 



the long distance physics in terms of collinear, soft and ultrasoft (usoft) fields. The short 
distance physics is absorbed into Wilson coefficients. So in SCET a quark field gets replaced 
by these different fields, which can each run around in loops. Currents in SCET often 
generate ^ divergences at one-loop that are removed by renormalization. We will see that 
such divergences show up in individual SCET anomaly diagrams, but cancel when they are 
summed. Also in SCET there are vertices with two quark fields and more than one gluon, 
which lead to non-triangle anomaly diagrams. These graphs turn out to be important. 

Which SCET fields one needs to consider, depends on the physical process one has in 
mind. We will be looking at SCETj with one collinear direction n'^. The fields in this theory 
are collinear quarks ^„ and gluons A'^ and usoft quarks qus and gluons A!^^. The momenta 
of the collinear fields scale like = {n ■ p,n ■ p,p±) ~ <5(A^, 1, A) and for usoft fields as 
Pus ~ QA^. Here Q is the hard scale, A the SCET expansion parameter and and are 
light-cone basis vectors satisfying = = and n-n = 2. One can for example take 
n'^ = (1,0,0,1) and n'^ = (1,0,0,-1). For a typical process A^ = Aqcd/Q- We will also 
need the power counting of the fields, which are listed in table [B 

To get a systematic expansion in the power counting, collinear fields have both a label 
momentum p containing the collinear momentum and a (usoft) coordinate x, ^n,p{x). The 
label momenta are picked out by the label operator P, e.g. V^^n^p = p'j_^n,p- The collinear 



^ In the Standard Model only the top would get integrated out before W and Z do. For simpHcity we 
assume the existence of a heavy electroweak doublet. 
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collinear 


usoft 


covariant derivatives 


operator 
power counting 


TT-An An^ Tl-An 

A AO A A2 


Q.US Aus 

A3 A2 


in-Dc iDc^ in-D iDus 
AO A A2 A2 



TABLE I: Power counting of the various fields and operators in SCET 



covariant derivatives are then defined as 

in-Dc = n-V + gn-An,q, iD^^ = + gA:^^g, in-D = in-d + gn-An,g + gn-Aus, C^) 
and the usoft covariant derivative as 



(2) 



Note that these covariant derivatives are each homogenous in the power counting. We are 
now ready to write down the leading order Lagrangian for a collinear quark 



Cfl =Ux)(in-D + i 



(3) 



where it is understood that the (suppressed) label momenta are summed over and that the 
label momenta of each term in L is conserved. 



This paper is organized as follows: we start in section UTAI by matching the chiral anomaly 
equation in QCD onto SCET, using the tree level relations between fields. For simplicity we 
first restrict ourselves to LO in the power counting, postponing the derivation up to NNLO 
to section IIII A[ Whenever we speak of e.g. LO and NLO in this paper, we will always 
be referring to the order in the power expansion A and not to radiative corrections. No 
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contributions are expected beyond one-loop due to the Adler and Bardeen theorem 

We then verify these SCET anomaly equations by computing the chiral anomaly from 
one-loop graphs in SCET. This is done at LO in sections FlI BIIII Dl NLO in section [Til Bl and 
NNLO in section HVl At LO only collinear fields play a role and so we expect agreement, 
since SCET with only collinear fields is just a boosted version of QCD. Things become more 
interesting at higher orders when diagrams with both collinear and usoft fields show up, but 
we find that the derived SCET anomaly equations are correct. We end with some concluding 
remarks in section IVl and list the full SCET anomaly equations up to NLO. The appendix 
contains a table of the loop integrals we used. 



II. THE SCET ANOMALY AT LO 



A. Matching the QCD anomaly onto SCET 



In this section we will calculate the chiral anomaly in SCET at LO. We start by expanding 
the left- and righthand side of the QCD anomaly equation 
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9' 



(4) 
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in the power counting, where J5 = ^^^■j^'if. We use the tree level relations between fields, 
effectively matching QCD onto SCET at tree level. In the next subsections we will explicitly 
check that the SCET anomaly equation we find, holds at one-loop. 

We introduce the following notation for the left- and righthand side of equation ([4]) 

J^d,4, J^^-J1- e^^.MF'-^Fn, (5) 

iOTT 

and write their expansions as 

J = JW + ^(5) + ^(6) + . . . ^ (6) 

where the order in A is indicated in brackets. When calculating the one loop anomaly 
diagrams, we will (as usual) restrict ourselves to two outgoing gluons. Gauge invariance 
then determines anomaly matrix elements for more gluons. Let p, q be the momenta of 
these two gluons. For simplicity we assume that they have no component in the _L direction, 
p± = q± = 0. We will also assume that the gluons are ±-polarized. The matrix element of 
is then given by 

{99\J\0) = ^^n■{p + q) {gg\n- MO) + '^zn-{p + q) {gg\n-J,\0). (7) 
The tree level relation between the QCD and SCET (collinear) field is 

A A2 

with the power counting as indicated. This leads to 

fl- ^ = Inilbin- (9) 

For n- J5 we get an t^. But t/i^n = 0, so the lowest non- vanishing contribution is at 0{\^) 

However, it is important to note that n-J^'^^ contributes at the same order as the n-J^ term 
in the anomaly equation because of the power counting for collinear momenta, 

{gg\J^'^\0) = ¥n-ip + q) {gg\n-4'^\0) + ^tn-ip + q) {gg\n-jt^\0). (11) 

We will first study these terms separately by imposing n-{p + q) = or n-(p + g) = in 
the following subsections. With these additional assumptions we have to keep gluons off- 
shell, 7^ 0, because everything would otherwise trivially vanish. At the end we will 
drop these additional assumptions, effectively combining the two cases. 
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Now we consider the righthand side of the anomaly equation. With our assumptions the 
important terms are 

in-D = in-d + gn-An^q + gn-Aus = in^+ ■ ■ ■ , (12) 
in-D = in-Dc + Win-DusW'^ = n^ + in^+ . . . , 

1 A2 

iD^^ = iDi^ + WiD^^W^ = gAi^^ + gA^^ + ... 

A A2 

W is the collinear Wilson line that only contains the n-A^ component of the gluon field 

W = l-g ^n-An + ... (13) 
n-V 

and can be dropped because we are only looking at ±-polarized gluons. Writing 

ea/^M.^"^^'^'^ = -j2 eap,.^D"^D''^DW, (14) 
the only contribution we need to consider at 0{X'^) is 

•^^'^ = -jf^ X Y Y tr[m-9 gA^>; n-V gA^,^^ - m-d gA^^^ n-V gA^^^ 

-n-V gA^^^ m-d gA^^ + n-V gA^^ m-d gA^^^ . (15) 
The corresponding matrix element is given by 

{gg\r'^'M = {n-p n-q - n-q n-p)5^^e^^^e^;[p)e^*{q), (16) 

where e^*(p), e^*(g) are the polarization vectors of the gluons and we assume the generators 
are normalized as tr [T'^T^] = . The two-dimensional e-*- is defined as 

i = |ea/3/.z.n°n^ = \i tr[75^7i7;^7^]. (17) 

Having derived the explicit form of the SCET anomaly equation at LO 

j(^) = (18) 

we will verify this result by calculating the one loop diagrams. 
B. LO with n-(p + g) = 

We start by considering the simplest case, namely (p + g) = for the momenta of the 
external gluons. The axial current is then given by 

{ggW'^'M = \in-{p + q) {gg\n-jf\Q) = \m-{p + q) {gg\lnil,U^) (19) 
and there are two diagrams that contribute at this order, shown in figure [H As mentioned 
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FIG. 1: Triangle and bubble diagram contributing to the anomaly at LO, for n-(j) + q) = 



before, taking the gluons on-shell would make everything vanish because oi n-{p + q) — 0. 
We therefore keep the gluons off-shell, which generally takes care of any IR divergences 
as well. Because we assume n- {p + q) = 0, we still need to be careful when combining 
denominators as (/ + p)^ and (/ — g)^. We will elaborate on this for the bubble diagram. 

In dimensional regularization with d — A — 2e, the triangle diagram is given by 



At^ \in-ip + q) (-1) Jdlti 



Iph , hip \ i .in-{l + p) 



n-l ^ n-{l-q) J 2 " 2 P 

+ (p,p,A) ^ (q,u,B), (20) 



n-{l+p) n-l J 2 2 {I + p)' 

where / dZ = / d'^l(27r)-'^. We will always suppress polarization vectors. To avoid any 
complications related to 75 in dim. reg., we never (anti)commute with 75 in our calculations. 
Since the external momenta have no _L-component, we can immediately replace — > 
^llg'f, where A = ci - 2 = 2 - 2s. We then use 



lllp = A, 7^^7,^7, = -(A - 2)7,^ = 2^7,^, 
lll^li = ^9up - (4 - A)7,^^^. 



(21) 



This leads to 



At = -ig'n-ip + qW^'e'^.ih + eh + eh - (1 + '2e)h) + (p, P, A) ^ (q, u, B) 

.2 



— n-{p + q)n-p 5 e^^. 



where 



n-(l — q) n-(l + p) !"]_ 
n-l {l-qf P {l+pf 
n-{l + p) l\ 
{l-qY P {l + pf' 



h- 
h 



dl 
dl 



n-{l-q) I] 



{l-qf P (1+pY' 

n-l ll 
{l-qf P {l + pf 



(22) 



(23) 



We computed these integrals using an appropriate Feynman parametrization. A table with 
all the necessary integrals can be found in the appendix. 
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The contribution from the bubble diagram (right panel of fig. [T|) is 



Ab = lin-{p + q) (-1) J m tr 



75 « 



•rjin-ij. — q) 
2 {I- qy 



X ig 



n-{l + p — q) 



+ 



ig'^n-{p + q)5 



AB 1. 



d/ — 



n-l n-{l + p — q)J {l — qY {l + pY 



t/l .r^n-{l + p) 
n-l I 2 ^2 {l+py 

n-{l — q) n-{l + p) 



(24) 



By sending / —l—p + q we see that the two terms are the same (giving a factor of 2, rather 
than canceling each other). We do not have a second contribution from (p, p, A) ^ (g, z/, B) 
for this diagram. As mentioned before, we need to be careful in dealing with IR divergences 
here. Writing 6 = n-{p + q) -^0, we encounter when doing this integral. Obviously 
the order of taking 5 — > and e — > matters and we should first take e ^ with 5 7^ to 
regulate any IR divergences. Proceeding along this way, we still find As = 0. 

We now add At and Ab and compare with (fi'fi'|^*^'*-'|0) as given in (fT6|l . They agree. This 
was expected since at this order only collinear fields are involved, which is just a boosted 
version of QCD. 



C. LO with n-(p + g) = 

We move on to the next case: n-{p + q) = 0. This time the relevant term of the axial 



current comes from n-J^'^^ 

(ggU'^'^lO) = ltn-{p + q) {gg\-^^ipi^^^^,^^ ipiUo). (25) 

in- D c in-JJc 

It can be part of a diagram in various different ways, for example a gluon can now come out 
of the current. These different ways are pictured in fig. [2] and the corresponding expressions 
are 

1) - -77^^75— ^ U 2) - e.-77^^75^P^ (26) 
n-[l + p) n-[l — q) n-[l + p) n-[l — q) 

n-[l + p) n-[l — q) n-[l+p) n-[l — q) 




FIG. 2: n-Jl^' can be inserted into diagrams in various ways 
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Up 



' P.P.A 



P, P, A 



C1) 



I 

q, V, B ^ - ^ ^ 

^ \ p, p,A 

(X) ^ nrnmp 



Up 




C2) 



X) 




P, P, A 



l-q 



D) 



q, V, B 




P, P, A 



FIG. 3: Diagrams contributing to the anomaly at LO for n-{p + q) =0 
First of all we get a triangle fig. [3]A and bubble diagram fig. [3]B coming from 1) 



Aa + Ab = tg'n-ip + g)5^^e;,((l + 2e)h - eh - eh -h~{l + 2e)h) 
+ (p,p,A) ^ (g, u,B), 



(27) 



where 



h 



m 



n-l {l-qf P (l + py 



n-{l-q) {l-qf P {I + pf 
I' 



, h 



/1 



n-{l + p) {l-qf P 

n-l Ij 

n-{l - q) n-{l + q) {I - q^ P (/ + pf 



(28) 



n-l {l-qY il + pY' 



We also have diagrams fig. [3lGl with 2) and fig. [31CI2 with 3), where one gluon comes out of 
the current 



Aci + Ac2 = 2ig^n-{p + q)5^''e%{h + eh) + (p, p, A) ^ (g, z/, 5), 



with 



d/- 



n-l 



(29) 



(30) 



'n-{l-q) n-{l + p) P {l+py ^' J ^" n-{l - q) P {I + pf 

Finally there is also a diagram fig. [3D where both gluons come out of the current, which 
involves 4). This turns out to be zero, Ad = 0. 

We simplify the remaining integrals using 

/2 = f -n-l n-l = {l- qf -n-{l - q) n-{l-q) = {I + pf - n-{l + p) n-{l + p), (31) 

where l']_ = — if g^" = diag(H ). For example 

II n-{l + p)ll 



d/ 



n-{l + p) {l-qY P {l-qfPil+py 



hA + hs- 



(32) 
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Doing the math, we find 

t 
57r 



{99\J^''^\^) = Aa + Ab + Aci + Ac2 + Ad = n-{p + q)n-p 5^^e^,. (33) 



Again this agrees with (5'5'|^^'*-'|0), as expected because the calculation so far only involved 
collinear fields. 



D. LO general case 

Finally we now drop the additional assumptions (e.g. n-{j) + q) = 0) and study the general 
case. These assumptions allowed us to look at one term of {gg\J'^^^\0) in (fTB at the time, 
so we basically have to add the anomalies of sections III Bl and III Ci The only other place 
these assumptions were used, is in simplifying the Feynman integrals. We can now take the 
gluons on shell and we will assume 

n-p = 0, n-q = 0. (34) 
After repeating the above analysis and performing the Feynman integrals we find 

{99\J^'^\0) = £^ n-p n-q 5^^e;t = (gglJ'^'^lO) . (35) 

Thus the SCET anomaly equation at LO J'^'^^ = J-'^'^\ which we derived by expanding, is 
correct at one-loop. In the previous sections the bubble diagrams in fig. [3]CI1, C2 contributed, 
but the bubble diagram with the two gluon vertex in fig. [U and fig. [3|B turned out to be 
zero. However, in this general case with on-shell momenta we do get a genuine contribution 
from them. 



III. THE SCET ANOMALY AT NLO 

A. Matching the QCD anomaly onto SCET at higher orders 

We will now move on to calculating the chiral anomaly in SCET beyond LO in the power 
counting, by taking into account the higher order terms in eq. (I6|). Again we start by match- 
ing the QCD anomaly equation onto SCET at tree level. The derivation is similar to that 
in section III Al but this time we need the tree level relation between the QCD and SCET 
fields to higher order in A. Using the result from we have 

tj/r,^ ^r^ + W^tprtw^^n + ... (36) 




in-Dr 2 " n-V 2 



A2 A3 



We find 



n- 


lm 




n- 


t(3) 


= 0, 


fi- 


t(4) 
'•^5 





(37) 
(38) 
(39) 
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and 

n-J^^ = il^-^ (40) 

in-Dc m-jJc 

n.jf) = -2^^ipi^^ l^iWq^s + zpT^ W^^r.) + h.c. (41) 

in- Dc ^ ■ ' ^ 

Expanding the other side of the QCD anomaly equation gives 

{gglJ^'-M = {n-p n-q - n-q n-p)5^^e;^,e7(p)e^*(g), (42) 



,2 



{gg\:F^'M = ^ n-p n-q 5^^ e^f'X (pYb i^) . (43) 



^^2 - f - ^ - -fj.u 
^2 



(gglJ'^'^lO) = {n-p n-q, - n-q n-pr)6^''e';,e^^*{p)e^^*{q), (44) 



where the subscript in pr denotes residual momentum, p!^ X^. We already checked at LO 
that the SCET anomaly equation found this way is correct, JT"*-^^ = J-'^^\ In the remainder 
of this paper we will do the NLO and NNLO calculations. 

When we go beyond LO, we obviously get contributions to the anomaly from higher 
order currents such as n-J^ . However there are also contributions from diagrams involving 
the leading current and insertion(s) of subleading Lagrangians. To be specific, the SCET 
anomaly equation at NLO and NNLO read 

T{j(^) ^/:«} + = T{^(^) ^/:«} + ^(5), (45) 

T{^(^) i +T{^(^) ^£(2) I ^^1^(5) ,£(1)} +^(6)^ ^46^ 

where Tjj'^^) iC^^^} = J d^x T J^^\0) iC'^^^x) etc. £(1) and £(2) ^efer to all terms in the 
quark and gluon action of the respective order. 

Let us consider the righthand side of these equations. Calculating the anomaly at one- 
loop corresponds to tree level diagrams for The only non-trivial diagrams come from 
insertions on external gluon lines. If one would include insertions on external gluon lines 
you would get contributions on both sides of these equations, which are equal because these 
insertions are not part of the loop. We will not consider such diagrams and so the right of 
(145]) and (SI reduce to J^^^) and J^^^). 



We list the subleading Lagrangians we need below [13|, llJ, Il5|, Il6|, Il7|, Il8|, ll 

4^ = i^nW)zfts-^ {WUpt^ U + i^mPtW)-^ ^^us{y^% in). (47) 

4? = (^"'^^)*--^ *f (^^^'^) + ^^-^-1 (W^^^^^n), (48) 

4? = ^"^^ ^^^W^^ns + h.c, (49) 
where ig^ = [in- Dc,ip^]. 
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B. NLO with + = 



At NLO things become more interesting, because both collinear and ultrasoft fields are 
involved. To get the right power counting one of the outgoing gluons must be collinear and 
the other usoft. The only additional assumption that makes sense here is 

n-{p + q) = 0, (50) 

because n-p ~ 1 ^ ~ n-q. 
The axial current gives rise to 

{gg\T{j^^^ + J(5)|0) = + {gg\ (T{n-jP +n-jf ) |0) (51) 

V in-Dc in-JJc 

( - 2a^r±^U 75(W^g.s + tpf^ W^^r^) + h.c)) |0). 

Said in words, we get contributions from diagrams with n-jf^, one C^^^ vertex and C'^^^ 
vertices, and diagrams with n ■ jf^ and only C'^^^ vertices. There is no term in for 
a _L-polarized usoft gluon, so the usoft gluon has to come out of the vertex or the 
n-Jg^-* current. This leads to the diagrams in figure [4] and [5l Obviously we do not have 
(p, p, A) ^ (g, z/, B) terms in this case, since one gluon is collinear and one usoft. As it turns 
out, we already computed all the necessary integrals in the LO n-(p + g) = calculation. 
Here we have the additional simplification that we can generally drop n-q for the usoft 
momentum g, because n ■ q <^ n ■ p,n ■ I. 

Au yimuu^r .o, ^ nmmr bi) b2) 



l + P ' P.P.A Uq U" q,v,B p, p, A 

FIG. 4: Diagrams from n-J^'^^ contributing to the NLO anomaly. The usoft gluon vertex is C^^'^ 



\ p, p, A , \ p, p, A 




q, V, B , , _ - 

l+p l+p q, V, B P, P, A 

(5) 

FIG. 5: Diagrams coming from n-J^ contributing to the NLO anomaly 

The various diagrams lead to contributions 

+ A^i^ = 2^g'n.p S^^ej^il + e)Ch - h), (52) 

+ All^^ = 2zg'n.p 6^^ej;,{l + e)hA, (53) 

A^r + A^c2'' = '^^a'n-p 5^^ei.(/6 + e/7), (54) 

A^^^ + Al!^^ = 0, (55) 
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where I2A was defined in (l32ll . Adding the pieces we get 

{gg\i:{j^^ iC^'^] + J^'M = A^i"" + ^^2'' + Ki"" + Am"" + 

f^NLO , aNLO , aNLO , aNLO 

= 1^ n-p n-q S^'^ej, = {g glJ'^'^ \0) . (56) 
So the SCET anomaly equation is correct at NLO too. 

IV. THE SCET ANOMALY AT NNLO 

We will pursue our calculation to one higher order in A. At this order the loop momentum 
can have either collinear or ultrasoft scaling. 

Let us start by observing that there is a freedom in what you call label and residual 
momentum. Consequently SCET should be invariant under (for example) the following 
reparametrization 

n-Pc ^ n-Pc + fi-k, n-pr ^ n-pr — n-k, (57) 

where pc is the (collinear) label momentum, pr the (usoft) residual momentum and k some 
constant usoft momentum. This ties together JF^"^) and J-'^^\ basically predicting the latter. 
It is easy to check that (i42ll + (l44l) satisfies this. We will also calculate the loop diagrams, to 
verify the SCET anomaly equation at this order. Since both gluons are collinear this time, 
either n- {p + q) = ov n- {p + q) = can be imposed. We will restrict ourselves to the 
former. As can be seen from the expression for JF^^) in (l44ll . we need to keep the residual 
momentum components of the external momenta 

P = Pc + Pr = {0,n-pc,pj = 0) + {n-pr,n-Pr,Pr = 0) (58) 

to get a non-vanishing expression. Our assumptions are 

n-{Pc + qc) = 0, n-{pr + qr) = 0. (59) 

The axial current gives 

{gglTij'^^^ \ +T{^(^) ^ ^|^(5) + ^(6)|o) 

= \in-{p + q) {gg\[i:^n-jf^ \ iC^^^ + Tjn- + 

= \in-{p + q) {gg\{T[Lii,in \ iC^^^ + T|e„^75a iC^'''^] + 

+ (en^75W^g«. + h.c.)^|0). (60) 

We cannot make a diagram using n-jf^ and only vertices, since the n-J^'^^ current has 
a collinear and an usoft quark coming out of it. Therefore the only contributions come from 
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\ 1(2) / 




P, P, -A 




P. P, -A 



/ + p 



' P, P. ^ 



C2; 



l-q 



Up 



P, P, >^ 



D) 



/-q 




/ + p 



P, P, -4 



/ + p 



U" q, V, S 



L"' p, p, A 



FIG. 6: Diagrams contributing to the anomaly at NNLO for n-{p + q) = 



There is a bit of a technical issue with loop integrals that we need to discuss here. At 
LO and NLO we had collinear loops for which we combined label and residual momenta 

I / d/. This was possible because these diagrams only involved n-Z^, Ic and 

n-lr. Now we can get diagrams that also depend on other components of the residual 
loop momenta e.g. l^. This requires us to do / d/^(...) separately, which yields zero 
for a collinear loop since all factors of are in the numerator. This is not the case 
when the loop momentum is usoft, because then l:^ appears in the denominator as 
well. Diagrams with collinear loops may still have non-zero contributions, for example 

J d/r fi-{lr + pr){-..) = n-pr^i^J d/r(---)- Bccausc wc take the external = = 0, 
many diagrams vanish immediately. 

First of all we get the same diagrams as at LO, but with an extra C^"^^ insertion (fig. 

A NNLO , A NNLO , a NNLO , a NNLO , a NNLO _ 

- 2ig^n-{p + q)n-pr 5^^e^, {h + 'ieh - (1 + 2e)i^ - h) + (p, P, A) ^ (g, z/, 5), (61) 



where 



h= I dl 



n-il + p) Ij 



n-l /I 



n-{l + p) {l~qf P (I+pY 



h 
h 



I' 



{l-qY P [i + pY' 

n-jl + p) l\ 
n-l {l-qf {1+pY 



(62) 



We also get a triangle (fig. [DCIl, C2) where one of the vertices is and a bubble diagram 
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(fig. EP) with a vertex 

aNNLO , aNNLO , aNNLO _ 
^Cl + ^C2 + — 

2tg^n-{p + q)n-pr S^^'ej^iefs - {I + 26)fg + {p, p, A) ^ {q,u,B) + ho), (63) 

with 

{l-qf /2 H+pY' ' ~ y n-(/ + p) {l-qf P [l + pf ' 
J {n-{l + p)Y 

= J ^'n-ni-qf ii+pr 

Finally there is also a mixed collinear-usoft diagram (fig. [6E) with an usoft loop momentum 



-n-q n-l n-p 
-n-q n-{l — q) P n-pn-{l+p) 



A^'^ = -^g'n.ip + q)S^^ej^ / d/ ^ J (66) 

(p, p, A) ^ (g, z/, B) 
= 0. 



In a naive calculation where you forget to treat / as a residual momentum that is multipole 
expanded in the collinear propagators, you would misleadingly find a non-zero result for 
this diagram. 



(67) 



Rewriting the loop integrals using (ISTll and 

d 1 n{l + p) d 1 n{l + p) 



d{n-p){l+pf {1+pY' d{n-p){l+pf {1+pY' 
we get 

i, = + n-p ^Jg, (68) 

h = h-Is + h- n-p ^/g. (69) 
This reduces the calculation to 

_ aNNLO , aNNLO , , aNNLO , aNNLO 

— ^41 "T ^49. -t- ■ ■ ■ -t- /In -|- /1r 
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= -2ig^n-{p + q)n-pr S^^'ej, (els - n-p + (1 + 2e)n-p - i/ 

+ ip,p,A)^iq,u,B). (70) 

Working out Jio, one finds zero. The remainder is fairly easy to evaluate if you take 
(p, p, A) (g, z/, i?) before performing the integral over Feynman parameters. We find 



^2 



= 1^ n-(p + g)n-p, ^^^e^ = (^?^?|^(^)|0). (71) 
Once again we see by direct computation that the SCET anomaly equation is correct. 
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V. CONCLUSIONS 



We derived the chiral anomaly equations in SCET up to NNLO by matching the QCD 
anomaly equation onto SCET, using the tree level relations between fields in QCD and 
SCET. Gathering all the expressions, the anomaly equations up to NLO read^ 



j(4) ^ jp(4)^ 



(72) 
(73) 



where 



J 



(4) 



in- D, 



+ 



d. 



in-Dr 



(74) 



J 



(5) 



+ 



d. 



en7^5(H^g.. + w-^ zpr^ w^^n) 



+ h.c., 



167r2 
1 

8^ 



e^j,tr[m-D in-D^ iD^^ iD^'^ ± permutations], 
e^^tr[m-D in-D^ iD^^ {WiD^^W^) ± permutations]. 



(75) 

(76) 
(77) 



Here the sum over permutations means the sum of all diff'erent orderings of the operators 
multiplied by the sign of the permutation. Note that J^'^^'^ has an additional factor of \ 
compared to !F^^\ which comes from the redundant interchange^ of the iDjr. The relevant 
terms of C^^^ can be found at the end of sec. IIII A[ 

By explicitly calculating one loop graphs we verified these equations. Although we have 
not checked every possible term in the operator relations, our work suggests their correctness. 
At LO the anomaly only involves collinear fields and SCET with only collinear fields is just 
a boosted version of QCD, so we expected agreement. It was not a priori clear how the 
anomaly equations would work out beyond that order. 

Of course there are still higher orders to consider, which might be worth pursuing if 
one has in mind applying the anomaly in some specific process. JF contains terms up to 
order A^, whereas J has terms of arbitrary high order. This can be seen from replacing 
l/{in-Dc) 1/ {in-Dc + Win-DusW^) in ( l36ll and expanding. If our result extends to higher 
orders, then certain matrix elements of the currents will have to cancel each other beyond 



^ We only studied some terms of the NNLO anomaly equations and would need the higher order terms in 

l(36l) to write down the full expression at NNLO. 
^ In our derivation of T'^^'^ in lfT6|) we did not have this factor ^ , because we were distinguishing the outgoing 

gluons. 
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APPENDIX A: USEFUL INTEGRALS 



All the loop integrals needed for our calculations are given below: 

S'l ,„ = ^^^i^ (1 + (log 4„-1)b + ...) + (m'f-- ( A2) 



/ 



n-{l + a) [P — m^y lovr^ r(z/) n-a 

"3 f — ^ r4^^ ^ = ^ (1 + (log 47r -!)£ + ...) — (A4 

n-{l + a) {P-m?Y IGtt^ ^ ^ ^ ' ' n-a ^ 

where / d'^/ = / d%27r)-'^ and d = 4 - 2e. Note that ll = -Z^ when gt"" = diag(H ) 
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